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Abstract 

The equations of hydrodynamics are rewritten in sense of function- 
als with values in Non- Archimedean field of Laurent series or Re- 
distributions. A new ideology for understanding of conservation laws 
is proposed. A set of nonlinear algebraic equations suitable for the 
numerical work is given. The Newton iteration method are used for 
calculation of microscopic shock profiles for the equations of compress- 
ible flow. 
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1 Historical remark and Introduction 

In 1943-44, von Neumann became convinced that the calculation of the 
flows of compressible fluids containing strong shocks could be accomplished 
only by numerical methods. He conceived the idea of capturing shocks, 
i.e., of ignoring the presence of a discontinuity. Employing a Lagrangian 
description of compressible flow, setting heat conduction and viscosity equal 
to zero, von Neumann replaced space and time derivatives by symmetric 
difference quotients. Calculations using this scheme were carried out; the 
approximation resulting from these calculations (see j|] ) showed oscillations 
on the mesh scale behind the shock. Von Neumann boldly conjectured 
that the oscillations in velocity represent the heat energy created by the 
irreversible action of the shock, and that as Ax and At tend to zero, the 
approximate solutions tend in the weak sense to the discontinuous solution 
of the equations of compressible flow. 

In it was counterconjectured that von Neumann was wrong in his 
surmise, i.e., that although the approximate solutions constructed by his 
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method do converge weakly, the weak limit fails to satisfy the law of con- 
servation of energy. 

In Q J.Goodman and P.Lax investigated von Neumann's algorithm ap- 
plied to the scalar equation 

u t + uu x = (1) 

(it is called the Hopf equation [Q]), in the semidiscrete case. Using numer- 
ical experimentation and analytical techniques the demonstrated the weak 
convergence of the oscillatory approximations, and that the weak limit fails 
to satisfy the scalar equation in question. 

Von Neumann's dream of capturing shocks was realized in his joint work 
with Richtmyer in 1950, see Q. Oscillations were eliminated by the judicious 
use of artificial viscosity; solutions constructed by this method converge 
uniformly except in a neighborhood of shocks, where they remain bounded 
and are spread out over a few mesh intervals. The limits appear to satisfy 
the conservation laws of compressible flow. The conservation of mass and 
momentum is the consequence of having approximated these equations by 
difference equations in conservation form; but the von Neumann- Richtmyer 
difference approximation to the energy equation is not in conservation form. 

In the paper T.Hou and P.Lax compared the results of a von Neumann- 
Richtmyer calculation with the weak limit of of calculations performed by 
von Neumann's original method. 

The aim of this paper, at first, it is to introduce the idea of understanding 
of conservation laws; the second is to proposed a numerical method for cal- 
culation of hydrodynamic shocks profile without using differences schemes. 

We believe that there is no discontinues in the nature of shocks. Dis- 
continues solution of hydrodynamic equations only a rough mathematical 
model of shocks. When viscosity is taken into account, for example, the 
shocks are seen to be smeared out, so that the pure mathematical surface 
of discontinuety are replaced by the thin layers (10 -7 -10~ 8 m) in which 
pressure, density, temperature and etc. vary rapidly but continuosly (see 
Fig§. 

The equations of compressible flow in one space dimension can be written 
in the following Lagrangian form: 

ut + Px = 

v t -u x = (2) 
et + pu x = 

It is a classical approach. Here u is velocity, p is pressure, v specific volume, 
and e internal energy, connected with p and v via an equation of state. Here 
we mostly study the case when the equation of state is given by the 7-law 
with 7 = 1.4 

pv 
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Substituting that into the (||) we get the following equation 

— ; (pv) t + pu x = 0. (3) 
7-1 

The first equation of (Q) is conservation of momentum, the second equation 
conservation of mass, the third the work equation. 




Figure 1: Steady-state plane shock 

We will consider a steady-state shock. Imagine a long pipe containing 
a fluid initially in equilibrium (thermally and mechanically), into which a 
piston is pushing from one end with constant speed, as shown in Fig.[l]. In 
the presence of dissipation the specific volume, v, and the fluid velocity, u 
are as shown by the curves. Our problem is to calculate the exact shape of 
this curves. The shock is steady, at least approximately, after it has gone 
to a sufficiently great distance from the initiating piston. Then u,v,e, etc. 
depend on x and t only through the combination y = x — ct, where c is a 
speed of the shock relative to the original, or Lagrangian, coordinates. 

Now, we are going to study the equations (||) in specific sense. Namely, 
we are going to rewrite mentioned equations in the sense of R(e)-distributions 
We give the definition of the special kind of solutions of the some conser- 
vation laws in the sense of R(e)-distributions and consider the method for 
the numerical calculations of the smooth shocks and soliton like solutions 
of the Hopf equation and equations of compressible flow in the mentioned 
sense. This method based on orthogonal system of the Hermite functions 
as a base for calculation of such solutions (i.e. shocks and infinitely narrow 
solitons). Calculations of profiles of infinitely narrow soliton and shock wave 
are reduced to the nonlinear system of algebraic equations in R n+1 , n > 1. 
We proved, using the Schauder fixed point theorem that the mentioned 
system has at least one solution in R n+1 . We showed that there is possibility 
to find out some of the solutions of this system using the Newton iteration 
method We considered examples and numerical tests. We also should 
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emphasis that proposed numerical approach do not use a difference scheme 
(see also ||). 

Let us consider a bit of theory which we will apply to conservation laws. 



2 Non- Archimedean field of Laurent series and Re- 
distributions. 

The theory of Non- Archimedean fields was considered in the book by A. H. Light- 
stone and A.Robinson 0. 

Definition 2.1 A Laurent series is a formal object 



EC c-™- 
C,n+k£ 



7 + k 

n=0 



where k is a fixed (i.e., fixed for this Laurent series), each £j £ R ; and either 
£fc 7^ or each £j = 0. 

The Laurent series ^^lo^™ 6 ™' wnere Co — 1 an d £ n = if n > 0, is 
denoted by 1. It is easy to see that the Laurent series is a field. Let us 
denote it by R(e). The norm on the field of Laurent series can define 

\x\ u = e _i/ ^ for each x 6 R(e) 

(in place of e can use any number greater than 1). The function v(x) is a 
Non-Archimedean valuation. Define 

(oo \ oo 

Zn +k e n+k j = k if Yl Zn +k e n+k + 0, i k + 0. 

The norm | • \ u have properties 

1. \x\ u = if and only if x = 0, 

2. — \%\v ' \y\ui 

3. \x + y\ v < m&x{\x\ u , \y\ u } . 

Here, we propose a general construction of the R(e)-valued generalized 
functions [1C], [M. These objects are a natural generalization of Sobolev- 
Schwartz distributions. We call them as R(e)-distributions. 

1. Consider all functions f(x,e) G C°°(R x (0, 1)) such that integrals 

o 

f(x, e)tp(x)dx 

exist for any e and for all ip{x) from a given class of functions X (X 
can be Cg°(R), §(R) and etc.). 
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2. Suppose also that J f(x,e)ip(x)dx is a number af j£ (rp) from the field 

— oo 

of Laurent series R(e). 

3. The two functions f(x,e) and g(x,e) call equivalent with respect to 
test functions X if and only if 

+00 +00 

J f(x,e)^(x)dx = a fiS (ip) = a g>e {^) = j g(x,e)^(x)dx. 

—00 —00 

The equality means in sense of the field of Laurent series R(e) for all 
functions tp € X. Classes of equivalent functions call R(e) -functions. 
The expression 

+ OC 

s ix 



J f(x,e)ip(x)da 



associates a number from R(e) with every ip. Such a quantity is called 
a functional. In this case a linear functional map X into the Non- 
Archimedean field R(e). Call these functionals as R(e) -distributions. 

Thus, 

Proposition 2.1 R(e) -function f(x,e) = if and only if 



i-oo 

//(„. W .)* = o 6 ae> 



for every -ip from X. 

The set of all R(e)-distributions denote by IR(X) 

Remark 2.1 Recall that the idea of representation of a function f 6 Lj 0C (R) 
in terms of a linear functional 

+00 



C£°(R) 3 V' — ► J f{x)ip{x)dx € R 



frased on well-known proposition that if f G Lj 0C (R) and / f(x)ip(x)dx = 

/or ant/ ip G C^°(R) £/ien / = almost everywhere. 

Let us consider an example of the R(e)-distribution. 
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Example 2.2 Take X = Cg°(R) and f(x,e) = (p(x/e), tp(x) £ C£°(R) 
then R(e) -distribution can write in the following form. 



+00 +00 +00 



tp(x/e)ip(x)dx = e J ip(x)dxip(Q) + e 2 J xcp(x)dx ^-^ - + . . . . 
—00 —00 

iVote f/mf (p{x/e) converges to the function 

{ ' ~\ 0, i/x^0. 

Last function almost everywhere equals to zero. 

Like Sobolev-Schwartz distributions we can differentiate R(e) -distributions. 
For example, 

+00 +00 

/d f d 

—ip(x/e)ip(x)dx = — ip(x/e)—ip(x)dx, 
dx J dx 



00 +00 +00 

ip(x/e)—tjj(x)dx = —e / (p(x)dxip' (0) — e 2 / x<p(x)dx- 



dx y ' J J 1! 

—00 —00 

It is evident that R(e)-distributions are more general objects than Sobolev- 



Schwartz distributions [14] 



3 Conservation laws. Non- Archimedean approach. 

A conservation law asserts that the rate of change of the total amount of 
substance contained in a fixed domain G is equal to the flux of that substance 
across the boundary of G. Denoting the density of that substance by u, and 
the flux by /, the conservation law is 

d 
dt 



u(t, x)dx = — I f ■ ndS. 
G JdG 



Applying the divergence theorem and taking d/dt under the integral sign we 
obtain 

/ (u t + div f)dx = 0. 
Jg 

Dividing by vol (G) and shrinking G to a point where all partial derivatives 
of u and / are continuous we obtain the differential conservation law 



ut(t, x) + div f(u(t, x)) = 0. 

Note, that if f(u) = u 2 /2 then we obtained the Hopf equation (||). In 
general, previous calculations based on the following well known proposition. 



6 



Proposition 3.1 If G £ L| 0C (R) and J G(x)ip(x)dx = for any ip £ 

— oo 

Co°(R) then G = almost everywhere. 

Definition 3.1 Let us consider two sets of the smooth functions, depending 
on a small parameter e £ (0, 1]. Let us take all functions v(t, x, e) which have 
the type 

fx — ct 

v(t, x, e) = Iq + Alip 



lo, Al, c are real numbers, Al ^ and (p £ §(R), J (p(y)dy = 1. We denote 

— oo 

this set of functions by I. We call I as a set of infinetely narrow solitons. 

Definition 3.2 Now, let us take all functions w(t,x,e) which have the type 

' x — at s 



w(t, x, e) = ho + AhH 
ho ,Ah,a are real numbers, Ah * and H(x) = J { y)dy, f { y)dy = 1 

— oo — oo 

and 6 £ §(R). We denote this set of functions by J. We call J as a set of 
shock waves. 

It is natural to consider conservation laws as an integral expressions 
which contain the time t as parameter. Therefore, we introduce the following 
concept. 

Definition 3.3 The function v £ I (or w £ J) will be a solution of the 
Hopf equation up to e~ l , I £ No in the sense of 'R(e) -distributions if for any 
t£[0,T] 



^ +oo 

/ {v t (t,x,e) +v(t,x,e)v x (t,x,e)}?p(x)dx = ^2^ k e k £~R(s}, (4) 

— rvn k = l 



+?° +0O 

/ {w t (t,x,e) +w(t,x,e)w x (t,x,e)}ip(x)dx = ^2r] k e k £TL(e) (5) 

-oo k=l 

for every ip £ §(R). In case when I is equal to +oo the function v(t,x,e) 
(or w(t,x,e)) exactly satisfies the Hopf equation in the sense of R(e)- 
distributions. 
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Certainly, one can consider instead of the Hopf equation some conserva- 
tion law. 

Prom mathematical point of view, we deal with a infinitely differentiable 



functions in definitions |3.1| and |3.2| , so that we avoid the problem of distri- 
bution multiplication. From physical point of view, functions from the set 
I or J can describe fast processes. Mathematical models of such processes 
based on functions from I or J may give additional information and take in 
account a short zone where physical system make a jump from one position 
to another. 

Thus, we will consider solutions of the Hopf equation which are infinitely 
narrow solitons or shock waves. It easy to see that 

v(t,x,e)^{ ! + AW0) ' % X Z Ct J as e-0 
i h, if ct. 

w(t,x,e) — ► ho + AhH(x — at), as e^O 
H is Heaviside function. 



4 Method for the numerical calculations of the mi- 
croscopic profiles of soliton like solutions of the 
Hopf equation in the sense of R(e)— distributions. 

Thus, conservation laws are integral expressions. Therefore, it is natural, 



that one can interpret the Hopf equation in the sense of the definition 3.3. 

We will seek a solution of the Hopf equation in the type of infinitely 
narrow soliton, i.e. let us v £ I. Substitute v(t,x,e) into integral expression 
(Q) using the following formulas 



oo 

d ( ( x — ct 



+oo -. 

Mx)dx = Y* ce k m k -^ k+1 Xct), (6) 



k=0 



oo 

8 { 1 2 ( x ~ °t 



dx\2 ip \ e 



\ «, +oo 

\U{x)dx = Y J -e k 9k^ {k+l \ct). (7) 
' > k=0 



We denote 

+oo +oo 2 

m k (<p)= J y k v(y)dy, g k (y)= j y k ^-d y , k = 0,1,2,.... (8) 

— 00 —00 
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Thus, we obtain 

+00 



r +£? j,(k+i) ( f \ 

/ {v t + w x } ^dx = { Al (c - k)m k - (Al) 2 g k ) e k W - 1 ' . (9) 

J k=Q 



From the last expression we have conditions for the function ip{x). Namely, 
9k(<p) ~ ^f-m k ((p) = 0, k = 0,1,2... . (10) 
From the first (k = 0) we have 



-00 



c -xr = - = U v 2 (*)dz- (id 

Al mo 2 J 

—00 

Hence, we can rewrite conditions (|l0|) as follows. 

+00 +00 +00 

ip 2 (x)dx- / x k (p(x)dx= / x k ip 2 (x)dx, k = 0,1,2.... (12) 



Now, let us prove the following lemma. 

Lemma 4.1 For any non-negative integer n exists such function ip € §(R), 
<p ^ which satisfies the following system of non-linear equations: 

+00 +00 +00 

x k ip(x)dx = \ x k ip 2 (x)dx/ / <p 2 {x)dx k = 0, 1, 2 . . . n. (13) 



Proof. First, we will seek function tp(x) in the following type: 

<p(x) = coh Q (x) + cihx(x) + . . . + c n h n (x), (14) 

where 

hk( x ) = —7= e ~ x 7/2 are Hermit functions. (15) 



V¥k\ 



Then we substitute the expression ( |i~4"| ) into conditions (^). After that 
we will have nonlinear system of n + 1 equations with n + 1 unknowns 
(co, c i, C2, • • • , c n ). We write this system by the following way. 



Ax = Ji[x), x = (cq,ci, . . . ,Cn) (16) 
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A is a matrix with elements 



-t*QO 

A kj = J x k hj(x)dx = (-iyi k y/2nhf\0), i = V-l, k, j = 0,1,2, ... n 

— oo 

N is nonlinear map such that 

^ n n 

= — - J2(N(k)x, x)e k = £ (17) 



\x\, 

1 11 fc=0 fc=0 



Vector efc = (eo, e%, . . . , e n ) such that e k = 1 and ej = for all j ^ k. N(k) 
are matrices with elements 



hoc 



Nij(k) = j x k hi(x)hj(x)dx, i, j, k = 0, 1, 2 . . . n (18) 

— oo 

and functions 



.r 



Note that functions f k (x) are continuous everywhere except x = and 
l/fcO?)! < 11-^(^)11 due to Cauchy-Bunyakovskii unequality. Matrix A is in- 
vertible for any n because of det(A) is a Wronskian for the linear independent 
system of Hermit functions ho(x), hi(x), ... h n (x) and 

det(A) = (2tt) -r W(h o (0), hi(0), . . . MO)). 



We can write the system ( |16D as 

n 

fk(x)A~ x e k = F{x) or x = ^ -1 (X(f)) = F(f). (19) 



x = 

fc=0 



Let us describe the function F : R n+1 i — ► R ra+1 . It is continuous except 
x = and bounded. Indeed, 

n n 

\\F{x)\\ < WA- 1 ]] ^ ||JV(fc)||, r n = WA-'W £ ||JV(fc)||. (20) 

fc=0 k=0 

Let us consider function 3sf(x). It is continuous function everywhere in R n+1 
except x = and, moreover, !N(R™ +1 \{0}) C LTi where II i = {z E R™ +1 : 
zo = 1} is a plane. Further A _1 (Ili) = H2 where II2 = {y 6 R n+1 : 

n 

aoji/j = 1} is another plane. 

k=0 



-OO 



A 0j - = y hj(x)dx = (-i) j V2TThj{0), i = v 7 ^!, j = 0, 1, 2, . . . n 



-00 
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Thus, we can consider the function F(x) which is defined on the convex 
compact set C n = Tl2C\B[Q,r n ] such that F : C n i — ► C n , where B[0,r n ] is 
a closed ball with radius r n . Function F is continuous on the C n because of 
C n . Now we can use J.Schauder theorem. 



Theorem 4.2 (Schauder fixed-point theorem [|13[]) Let C be a com- 
pact convex subset of a normed space E. Then each continuous map F : 
C i — ► C has at least one fixed point. 



Hence, we can conclude that our system fll9|) and therefore system (jig ) 
has at least one solution. Thus, there is a function <p(x) which satisfy to 
conditions fll~3|) proposed lemma. 



Remark 4.1 Let us a function <p(x) satisfies lemma condition. If (5 G R 
then the function ip(x + (3) also satisfies lemma condition. Moreover, if 



oo 

2 



if (x) dx = a then if (ax) satisfies lemma condition. 

Thus, we can formulate the following result. 
Theorem 4.3 For any integer I there is a infinitely narrow soliton type 



solution of the Hopf equation (in the sense of the definition 3.S ) up to e 
with respect to the norm | • \ v , i.e. 



v(t, x, e) = Iq + Alip ( 



V £ 



-t-o© 

l ,Al,c are real numbers, and V e S(R), / ^y)dy = 1. Moreover, 



C-Iq 1 



Al 2 



tff(x)dx. (22) 



For example, calculations in case I = 7 give the "profile" f(x) (see Fig. 
^) for the infinitely narrow soliton v(t,x,s) = ip (^-^): 

= U c 2 (4^-2) c 4 (16x^ -48^ + 12)1 ^ 

where c = 0.66583, c 2 = -0.23404, c 4 = 0.05028, c = 0.25032 (c is a 
velocity of the soliton). Numbers Co, c 2 , c 4 and c were found approximately 
by iteration method using the following sequence. 

x m+ i = A- 1 (N(x t m )), m = 0,1,2,... . (24) 
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-4 -3 -2 -1 1 2 3 4 -4 -3 -2 -1 1 2 3 4 



Figure 2: The case I = 7, c = 0.25032. Figure 3: The case I = 13, c = 0.35442. 




-4 -3 -2 -1 1 2 3 4 -4 -3 -2 -1 1 2 3 4 



Figure 4: The case i = 15, c = 0.38267. Figure 5: The case I = 17, c = 0.40892. 



Matrix ^4 and a nonlinear N were introduced in the lemma proof. 

Calculations of soliton-like profiles f(x) for the Hopf equation in case 
I = 13, 15, 17, 19, 21 give us pictures (Fig. g, £§ |, §, g). 

For the Z greater than 21 matrix A is close to singular and calculations 
can be inaccurate. 



5 Calculations of the microscopic profiles of the 
shock wave solutions of the Hopf equation in 
the sense of R(e}— distributions. 

A solution of the Hopf equation in this case we will seek in the set J. Namely, 

' x — at" 



w(t, x, e) = ho + AhK 
ho, Ah, a are real numbers, Ah ^ and 



£ 



K(x) = J e{y)dy, J e(y)dy = 1, 9 G §(R). 
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0.8 
0.6 
0.4 
0.2 

-0.2 





-4 -3 -2 



-4 -3 -2 -1 



Figure 6: The case I = 19, c = 0.43357. Figure 7: The case Z = 21, c = 0.45678. 

Substitute x,e) into the integral expression (|5|) using the following 
formulas 



k=0 



Of 



k\ ' 



(25) 



R x-at \ d \ fx-at 



e dx 



N +00 

> tp{x)dx = £ 

' k=0 



rk- 



k\ 



(26) 



We denote by 

+oo 

>n k (9) = [ y k 6(y)du. ,-,(())- / x k 0(x)\ j <H, ,)<)<, \ <ix. /,•-(). 1. 2, .. . 

(27) 



Thus, we get 



/ T +oo 
{w t + ww x }4>dx = ^2{(Ah) 2 r k - Ah(a - h )m k } e 
k=o 



ki^!M 



From the last expression we have conditions for the function 9{x) 

r k (6) - ^^m k (6) = 0, fc = 0, 1,2... . 
An 

From the first (k = 0) we have 

a -^= f 9(x)l f 6{y)dy\dx = \. 



(28) 



(29) 



(30) 
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Therefore, we can rewrite (^) as 

+00 +00 



x k V[.c 



)dx = j x k 9(x) I J 6(y)dy J dx k = 0,1,2... . (31) 



The same method one can prove that there is such function 6{x) G S(R-) 
which satisfies the following conditions 

+00 +00 
i j x k 6(x)dx= j x k 6(x) \ j f)\ (/)(]// j d.r k - (). 1.2 . . . //. \:12 } 



Thus, we can formulate next result. 

Theorem 5.1 For any integer I there is a shock wave type solution of the 



Hopf equation ( in the sense of the definition 3.S) up to e 1 with respect to 
the norm \-\ u - 

w(t,x,e) =ho + AhK ( J , (33) 



ho, Ah, a are real numbers, Ah 7^ and K(x) = j 6{y)dy, j 9(y)dy = 1 

—00 —00 

and 9 £ §(R). Moreover, 

a — ho 1 



Ah 



(34) 



Note that the condition fl34j) is Rankine — Hugoniot condition for the 
velocity of a shock wave. 

As in previous section we seek function 0(x) in the following type: 

<p(x) = a h (x) + a\hi(x) + . . . + a n h n (x), (35) 



where hk{x) are Hermite functions. Calculations in case / = 7 give the fol- 
lowing "profile" ( 
h = 0,Ah = 1). 



x — at 

lowing "profile" (K(x)) for the shock wave w(t,x,e) = K ( ) (where 



K(x) 



j c 2 (4V - 2) ^ >-4** + 12) 1 ^ 



where c = 0.79617, c 2 = -0.53004, c 4 = 0.17923, c = 1/2 is a velocity of 
the shock wave (see Fig. ||). Numbers cq, c 2 , C4 were found approximately. 
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Note that the function K(x) is not unique. There is a different function 
K\(x) which satisfies mentioned above conditions. It has the following type 



KAx) 



+ 



CO 



+ 



ci2r 



+ c^-2)} e - rV2rfT+ 



c 3 (8r 3 - 12r) c 4 (16r 4 - 48r 2 + 12) i 



(37) 



^233!^F 



7T 



where c = 0.18357, c x = -0.73567, c 2 = 0.74733, c 3 = 0.15327 c 4 = 
—0.29539, c = 1/2 is a velocity of the shock wave (see Fig. ||). Coefficients 
Co, ci 3 C2, C3, C4 were found approximately by the Newton iteration method. 
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Figure 8: Graph of the function K(a 



Figure 9: Graph of the function K^x). 



Taking in account the Rankine — Hugoniot condition (p4[) we also have 



graphs (Fig. |10| |ll|) shock profiles. 
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Figure 10: First shock profile 1 - K(x). Figure 11: Second shock profile l-Ki(x). 

Here we describe how it is possible to find coefficients cq, c±, . . . , c n in this 
case by the Newton iteration method for the following system of nonlinear 
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equations. 



P(c) = Ac - 2^2(S(k)c, c)e k = 0, c = (c ,ci, ... ,c n ) 



(38) 



k=0 



Vector effc = (eo, ei, • • • , e n ) such that e k = 1 and = for all j 7^ fc. S'(fc) 
are matrices with elements 



Sij(k) = / x k hi(x) ■ / hj(y)dydx, i, j,k = 0,1,2 ... n (39) 



Matrix A have elements 

+00 

A 



00 



x l hj{x)dx, i, j = 0,1,2 ... n 



We can write the formula for the Newton iteration method 



(40) 



x m+ i =x m - [P'(x m )] 1 [P{x m )\ , (41) 
where [P'(x)] is a linear map depending on the vector x. 



[P'{x)] [h]=Ah-2\ Y^{S{k)x, h)e k + ^{S T {k)x, h)e k 



,k=0 



k=0 



(42) 



Calculations of shock profiles K(x) for the Hopf equation in case I = 8, 
9, 10, 11, 12, 13 give us the following pictures (Fig. [l|, |T|, ||, ||, ||, 0). 
Here, we show only two different types of the shock type solutions of the 
Hopf equation. We can find more solutions if we take a different initial data 
for the Newton iteration method. 
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Figure 12: Shock profiles when I = 8. Figure 13: Shock profiles when I = 9. 
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Figure 14: Shock profiles when / = 10. Figure 15: Shock profiles when I = 11. 
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Figure 16: Shock profiles when I = 12. Figure 17: Shock profiles when I = 13. 



Remark 5.1 It is not easy to see that there is exist function 

oo 

0{x) = ^2 a n h n (x), a = (a , oi, . . . , a n , . . . ) G Z 2 , 



(43) 



n=l 



such that 

+oo 



^ J x k 6{x)dx= j x k 6{x) I J 6{y)dy J dx, k = 0,1,2... . 



(44) 



We think that it is true. 



6 Calculations of the microscopic profiles of the 
shock wave solutions of the equations of com- 
pressible flow in the sense of R(e)— distributions. 

Now, we are going to study the equations of compressible flow 

u t + Px = (45) 
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v t -u x = (46) 



! (pv) t +pu x = (47) 



7-1 

in specific sense. Namely, we rewrite mentioned equations in the sense of 
R(e)-distributions. 

We will seek for a solution of the equations in the following form 

u(t, x, e) = u + Aull ( - — —J, (48) 



no, Am, v are real numbers, An / and U(x) = J U(y)dy, J U{y)dy = 1 

— oo — oo 

and U e S(R). 

p(t,x,e) = p Q + A P p(?—^j , (49) 

x +oo 

Po, Ap,a are real numbers, Ap / and P(x) = J P(y)dy, J P(y)dy = 1 

— oo — oo 

and P e S(R). 

fx — at \ 

v(t, x, e) =v + AvV i — - — J , (50) 

x +oo 

vq, Av,a are real numbers, Av / and V(x) = J V{y)dy, J V{y)dy = 1 

— oo — oo 

and V G §(R). Note that a is a velocity of the shock waves. 
In the other hand, we suppose 

U(x) = a h (x) + aihi(x) + . . . + a n h n (x), a = (a , a±, . . . , a n ), (51) 
P(x) = b h (x) + bihiix) + . . . + b n h n (x), b = (b , 6i, . . . , b n ), (52) 

V{x) = coh (x) + c\h\{x) + . . . + Cnh n (x), c = (cq, a, . . . , Cn), (53) 

where h^{x) are Hermite functions. 

We understand the solution of the system in sense of R(e)-distributions. 
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Definition 6.1 Functions u G J , p G J and v £ J is a solution of the 
system (45), ftjdj), fflj) up to e~ l , I G Nq in the sense of R(e) -distributions 



if for any t G [0,T] 



+?° +00 

/ N + Px} V(x)(ic = ^ i k e k G R(e) , (54) 



k=l 



/+00 
{«t - u x } ip(x)dx = (k£ k G R(e) (55) 



y j~7y {ptv + pvt} +Ptu x \ i>{x)dx = ^^rj k e k G R(e) (56) 

-co k=l 

for every ip G S(R). 

In case when I is equal to +00 functions u(t, x, e), p(t, x, e) and v(t, x, e) 
exactly satisfies the system jjU^), (ffij), ( fflO ^ n the sense of 'R(e) -distributions. 

Substituting u, p and v into ( p^) , (|55|) we get the following relations for 
the moments. 

-aAum k (U) + Apm k (P) = 0, fc = 0, 1, 2, . . . ra, (57) 



-aAvm k (V) - Aum k (U) = 0, fc = 0, 1, 2, . . . n. (58) 
We denote as usual by 

+00 / x \ 



m k {U) = / x k U(x)dx, m k (UU) = / x k U(x) / U(y)dy dx, (59) 



where k = 0, 1, 2, . . . . 

It is easy to find a from ( |57|) and (|57| ) when k=0. Indeed, 



A P Au fart 

a = ^ — > a = — —— (60) 
Au Av v ; 

These are Rankine — Hugoniot conditions. Indeed, from the last expressions 
we have 

- 2 = (61) 
At; 

See J. von Neumann and R.D. Richtmyer || formula (21). We also conclude 
that An < 0, Ap < and Av > 0. 
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Because of Am, Ap and Av some real numbers, therefore, all three vectors 
with coordinates mk(U), rrik(P) and mk(V), k = 0, 1,2, ...n, respectively 
should be collinear. However, 

mo(U) = mo(P) = mo(P) = 1. 

Hence, a = b = c. 

Substituting u, p and v into ( |56| ) and taking in account the last equalities 
we get the following relations for the moments (k = 0, 1, ... n) : 

a(Apv + Avpo) \ ~ J 2aApAt;\ ~ 
PoAn ^— 1 mfc((7) + j ApA-u — - j m k (UU) = 0. 

(62) 

When fc = we will have Rankine — Hugoniot conditions for our system 



{p A. - '(WAV*)} + 1 | ApA , _ ^ApA.j = Q 



(63) 



Using the equalities (|60|), (|63[), we will have 

, A P . 1 / . ,x A n Pi (7 + 1H ~ (7 ~ 1)^1 
7 • + -T- • v + - • (7 + 1) • Ap = or — = - — — — . 

Av 2 po (7 + l)vi - (7 - l)i>o 

(64) 

Then taking in account the condition (|6T|), we get 



Hence, 



I wo + 5 • (7 + 1)A«| a 2 = 7p . (65) 



V2 , 1/2 

7-^ • (66) 



(7 + l)-^ + 2y V vo 



Thus, we get the same formula for the shock velocity as in the paper of J. 
von Neumann and R.D. Richtmyer || (formula (72)). 
Thus, from (62) follows that 

m k (U) = 2m k (UU), k = 0,1,2, ... n. 

This system we already know how to solve by the Newton iteration method. 
See conditions (|32| ) and solutions in this case. Thus, we can formulate the 
following result. 

Theorem 6.1 For any integer I there is a solution of the system of equa- 
tions $4$), Wk ) i n the sense of the definition 6A if 

Ap Au 

a = "a ' a = "a - ' 
An Av 

jpoAu - — ?-j(Apvo + Av Po )^ + ~ |ApAn - -^-ApAvj = 0. 
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Figure 18: Velocity and volume. 



Figure 19: Shock profile of pressure. 



If we assume that vq and Av are known then we can calculate all con- 
stants Ait, Ap, a and profiles of the shocks. Let us take vq = 0.2, Av = 0.8. 
We denote by v\ = vq + Av, v\ is a volume before the shock formation. 
The quantity r\ = vi/vq is a measure of the shock strength. Denote by 
Pi = Po + Ap, v\ is a normal pressure (it is known, p\ = 1 atm) before the 
shock formation. From the ( |64"D we can get 

27£>i 



Ap 



7 — 2vq/ Av — 1 
Therefore, we can calculate step by step 



A P A 
-- — , Au = —a ■ Av. 

Av 



Shock profiles of the considered system (|45[) , (|46|) , ( |47|) one can find on 
pictures (Fig. 18, We considered only one type profile of the shock, the 
case when I = 13, vq = 0.2, Av = 0.8, p\ = 1 and 7 = 1.4. 

Finally, Ap = -28pi, a = 5.9, An = -4.7, a = (0.47894, 0.70727, 
0.31783, -0.4706, -0.53605, 0.20792, 0.4099, -0.057108, -0.16832, 
0.0088166, 0.029395). 

In conclusion we should emphasis that our calculation method looks 
like the Fourier method for linear differential equations but applied to the 
nonlinear equations. Our method allowed to obtain all known formulas for 
the shocks characteristics and, in addition, find a microscopic behaviour of 
shocks in the thin layer. According to the our model the pressure in the thin 
layer can be negative. It is possible that the concept of "pressure" in the 
thin layer (where the jump of p(t, x, e) took place) one should understand 
in special sense. We think that the role of e in the calculations can play so- 
called the average length of free movement of gas molecules. From our point 
of view the phenomena of the formation and development of characteristics 
of the shock wave should describe in terms of non-Archemedean distance or 
non-Archemedean geometry. 
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In addition, we can use Laguerre functions or harmonic functions in our 
calculations instead of Hermite functions. 

We hope that our approach will be useful for the problems of nonlinear 
Optics and problems of Quantum Mechanics. 

Research is partially supported by Belarussian Fundamental Research 
Foundation Grant No F99M-082. 
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